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INTRODUCTION 


Several operational aircraft have accumulated significant 
flight time at speeds sufficient to produce severe aerodynamic 
heating (refs. 1 to 3) . Even with this experience, there exists 
a lack of understanding of how accurately thermal stresses can 
be predicted on a complex structure. The ability to predict ther- 
mal stresses accurately has great impact on both the safe magni- 
tude of stresses and the long-term effect of thermal cycling on 
the structure. An effort to study how well thermal stresses 
could be predicted in detail was begun in reference 4 and con- 
tinued in reference 5. This paper provides the third and last 
part of this particular study about how laboratory measured ther- 
mal stresses compare with thermal stresses predicted using NASA 
structural analysis (NASTRAN) computer models. 

A test structure representing a portion of a hypersonic 
vehicle (ref. 6) was instrumented with strain gages and thermo- 
couples. This test structure was then subjected to laboratory 
heating representative of hypersonic flight conditions. Two 
finite element computer models of this structure were developed 
using bar and shear panel elements of the NASTRAN program (ref. 7) . 
Temperature inputs from the tests were used to predict model 
thermal stresses, and these were correlated with the test meas- 
urements . 



DESCRIPTION OF TEST SPECIMEN AND INSTRUMENTATION 


A test specimen resulted from work directed toward studying 
the feasibility of a hypersonic research airplane (ref. 6) . One 
structural concept resulting from this study was a heat sink type 
of structure. The basic philosophy was to use a variable- 
thickness, load-carrying skin to absorb heat as required to main- 
tain a certain skin temperature. The skin was supported by a 
titanium zee substructure as shown in figure 1. A beryllium- 
aluminum composite metal (ref. 8) was selected as the skin 
because of the light weight and large thermal capacity of 
beryllium. The substructure frames were a zee-type configuration 
with titanium being used as the zee material. A photograph of 
the substructure frames is shown in figure 2. A photograph of 
the complete specimen is shown in figure 3 where the skins are 
attached to the substructure. 

The location of the temperature and strain instrumentation 
is shown in figure 4. Seventeen thermocouples and seventeen 
strain gages were used for analysis purposes in this paper. Five 
thermocouples and strain gages were located on the skins, while 
the remaining twelve were located at strategic points on the 
substructure. Thermocouples and foil- type strain gages were 
used in the test temperature environment, which ranged from 
232° K (-42° F) to 633° K (679° F) . 


TEST PROCEDURE 


The general test setup is shown in the schematic in figure 
5. The specimen is completely encased but is supported so that 
it is unrestrained — in other words, the edges are free to rotate 
and translate in the plane of the skin. A rack of radiant heaters 
(ref. 9) is located such that the skins can be heated on the side 
away from the frames. A blower system is situated such that a 
mixture of air and gaseous nitrogen can be blown over the specimen 
if cooling is desired. 

A flight profile for a hypersonic mission is shown in 
figure 6. The profile is of a 6-minute rocket-powered mission 
(ref. 6) with a 1-minute cruise at Mach 6. Skin temperatures 
were calculated based on the computer program identified in 
reference 10. 

A time history of skin temperature of the profile was used 
to perform the test. The test specimen is cooled prior to 
heating to simulate the cold soak condition that occurs when 
rocket powered aircraft are air launched at high altitude. 
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After the specimen is cooled, the heating is started by controll- 
ing the skin temperatures with radiant heaters according to the 
prescribed temperature time history. 


DESCRIPTION OF THE STRUCTURAL MODELS 


An initial NASTRAN model (shown in figure 7) was developed 
for one of the symmetrical quarters of the structure. This 
model, designated ASTEEL, was a 97 element model constructed of 
NASTRAN bar and shear panel elements. The 97 element ASTEEL 
model was composed of 73 bar elements and 24 shear panel elements. 
In constructing the model, 50 grid points were used with 231 
degrees of freedom, 69 single point constraints, and 250 bulk 
data cards. Since the ASTEEL model was constructed of very long 
elements, an additional model, designated CSTEEL, was developed 
with a significantly changed grid arrangement. The CSTEEL model 
(shown in figure 8) contained 129 bar elements and 54 shear panel 
elements. In constructing the model, 76 grid points were used 
with 387 degrees of freedom, 69 single point constraints, ancf 
402 bulk data cards. The NASTRAN models, which represented one 
symmetrical quarter of the test specimen, each have a surface 
area of 0.34 square meter (3.68 square feet) and a nominal depth 
of 0.10 meter (4 inches). The grid numbering system common to 
both the ASTEEL and CSTEEL models is shown in figure 9. 


RESULTS 


The results of the hypersonic heating simulation will be 
examined from the viewpoint of the thermal stress distribution 
at a single instant of time and from the viewpoint of the time- 
history of thermal stress at discrete points. 

Time histories of temperatures and thermal stress measured 
during the laboratory heating test for each of the locations 
shown in figure 4 are presented in figure 10 as the solid lines. 

The time histories of temperature generally begin at a low value 
representing a high altitude soak condition. The temperatures 
then rise to some higher value reflecting a maximum Mach 6 heating 
environment. The time histories of thermal stresses are compared 
with the ASTEEL and CSTEEL NASTRAN models. The variation in 
the calculated thermal stresses results from longitudinal tempera- 
ture variations (Y-direction) . Since the ASTEEL model has only one 
element in the Y-direction, there is no way of accounting for 
temperature variations in that direction. The CSTEEL model has 
three elements in the Y-direction, hence, temperature variations in 
the Y-direction could be considered. 
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Thermal stresses in the skin are presented in figures 10(a) 
through 10(e). Frame thermal stresses are presented in figures 10(f) 
through 10 (q) . It can be seen that the correlation of the 
measured test data with the NASTRAN models is highly varied. 

There is significant difference between the measured and 
calculated thermal stresses at locations such as is shown in 
figures 10(f) and 10(1). There is very close agreement between 
the measured and calculated thermal stresses in other cases 
such as shown in figures 10(h), 10 (m) , and 10 (n) . 

The distribution of measured thermal stress is shown in 
figure 11 for four different instants of time. The distribution 
in the left frame and the right frame is shown in the upper plots. 

The distribution in the skin is shown in the lower plot. The 
measured thermal stress data is compared to both the calculations 
from the ASTEEL and CSTEEL models. 

Temperature Smearing 

A general problem is becoming apparent in the analysis of 
large structures for thermal stress. The amount of temperature 
information seems likely to lag the structural input requirement. 

The test structure in this paper is a prime example of this 
problem. The ASTEEL model has 50 grid points and the CSTEEL model 
has 76 grid points. Each of these grid points requires a temp- 
erature input. There were available only 17 thermocouples at the 
strain gage locations plus two other control thermocouples (see 
figure 12) . The available test data supplied only a small per- 
centage of the temperature information potentially required as 
direct input to the grid points of the structural model. Somehow, 
the 19 temperatures may have to be expanded to 50 in the case of 
the ASTEEL model and to 76 for the CSTEEL model. The technique for 
deducing the temperatures at the grid points for which no direct 
temperatures are available is known as temperature smearing. 

Temperature smearing is necessary when the temperature 
distribution of the structure is non-uniform, unsymmetrical , or 
otherwise incomplete. The temperatures for the test structure 
in this paper were significantly non-uniform as can be seen in 
figure 13. The numbers in the panel squares represent the temp- 
erature at the center of each of the panels at the four minute 
time slice (see figure 6). It can be seen that no symmetry exists 
in the temperature distribution of the specimen. Hence, a 
problem exists when the model represents only one quarter of a 
structure that does not have an identical (or mirror image) 
temperature pattern in the other three quarters. 

The typical manner in which temperatures were smeared is 
shown in figures 14 and 15. In figure 14(a), the temperature 
along the longitudinal line connecting grid points 25, 125, and 225 
and 27 are shown. Along this line there were only two temperature 
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measurements recorded during the laboratory test; at 
grid points 25 and 225. These measured temperatures are 
represented by the two circular symbols. The temperatures for 
the ASTEEL model were considered to not vary in the longitudinal 
direction (from grid point 25 to grid point 27), hence, the 
temperature distribution for the ASTEEL model is represented 
by the straight, dashed line. The solid line passing through 
the two measured points represents the temperature smearing used 
for the CSTEEL model. A similar plot is shown in figure 14(b) 
for a line connnecting grid points 1, 101, 201, and 3. Along this 
line there were only two temperature measurements; at grid points 
1 and 201. The dashed line represents the temperature used for 
the ASTEEL model. The solid line represents the smeared temp- 
erature distribution used for the CSTEEL model. 

The lateral temperature smearing is shown in figure 15. 

It can be seen that several measured temperatures are available 
along the line connecting grid points 37 and 1. The distribution 
of temperature input to the ASTEEL and CSTEEL models is identical 
along this line. There is much less temperature data along the 
line connecting grid points 237 and 201. The CSTEEL model is 
smeared through the two available temperature data points. The 
ASTEEL temperature distribution (dashed line) is seen to vary 
significantly from the distribution used for the CSTEEL model. 

This is because the ASTEEL temperatures along grid points 237 
to 201 are a direct extension of the temperatures along grid 
points 37 to 1. 

More sophisticated smearing techniques can be established. 

The basic measured temperature data can be smeared to other grid 
points using supporting calculative techniques. However, for this 
test specimen, there was little justification for more effort 
than the linear interpolations and curve fitting techniques that 
were used. 


DISCUSSION 


The primary purpose of this paper has been to evaluate how 
well thermal stress could be predicted on a heated laboratory 
test structure constructed of dissimilar materials. Predictions 
were generated using the finite element NASTRAN models ASTEEL 
and CSTEEL. The manner in which temperatures were smeared was 
found to have great impact on the correlation of test data with 
predictions . 

The basic correlation data is presented in figures 11(a) 
through 11(d) in a generally condensed form. The stresses in 
the two spars were predicted quite well except near the interface 
with the skins. The measured stresses near the skin juncture 
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on the right frame are near 210 MN/m (30 Klb/in 1 while the- 
NASTRAN calculated stresses are less than 35 MN/m 2 (5 Klb/in 2 ) . 

A similar situation exists for the left frame, however, the 
discrepancy is not as large in magnitude. The two structural 
models (ASTEEL and CSTEEL) produced very similar answers in the 
frames for all time segments. 

Thermal stresses in the skin compare closely with the 
NASTRAN models for the zero and two minute time segments. At 
time = 4 minutes the thermal stresses calculated with the two 
finite element models do not agree with each other and they do 
not agree with the measured thermal stresses.- It can be seen 
that variations as high as 62 MN/m 2 (9 Klb/in 2 ) between the 
ASTEEL and CSTEEL models exist. This reflects the problem 
associated with smearing the temperature adequately. It can ' 
also be seen that the thermal stress data calculated using the 
NASTRAN models does not agree closely with the measured thermal 
stresses at time = 4 minutes. The calculated thermal stresses 
at time = 6 minutes also disagree somewhat with the measured 
data, however, there is closer agreement between the two models. 

It is apparent from the results presented in this paper and 
the results of reference 5 that the thermal stresses calculated 
using the NASTRAN models are very sensitive to minor variations 
in the temperatures input to the node points. It was shown that 
smearing temperatures in more detail, as was done for the CSTEEL 
model, does provide a different answer from the more coarse 
temperature smearing of the ASTEEL model. Insufficient informa- 
tion is available to determine the extent to which additional 
temperature information would aid in comparing with the measured 
values. The general comparison of calculated and measured 
thermal stresses is quite good and there is little doubt that the 
finite element approach provided by NASTRAN results in correct 
thermal stress calculations. There were some local discrepancies 
that could be attributed to either inadequate temperature infor- 
mation or insufficient NASTRAN elements in the local area. 

The areas in which the NASTRAN models failed to predict the 
thermal stresses accurately were in the frame-skin juncture areas 
and in the skins themselves. These inadequacies could have great 
impact on the long term life of an airplane. The large measured 
thermal stresses at the frame-skin juncture are most pronounced 
at time = 0 minutes. These stresses most likely result from the 
contraction of dissimilar materials (i.e. the titanium frames and 
the beryllium/aluminum skins) at this low temperature soak condition. 
The general inability to predict thermal stresses at this point 
is apparent all through the time-history of the hypersonic test 
(see figures 10(f), 10(1), and 11). Large and unpredicted thermal 
stresses in this particular area could have great impact on the 
ultimate life of fasteners used at this juncture. Local buckling, 
yielding, and fatigue problems could cause serious problems in 
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the fastener area. It is felt that the discrepancies occurring 
at the skin/frame juncture were most likely due to insufficient 
modeling elements rather than temperature problems. 

The discrepancy between measured and calculated thermal 
stresses in the skin area was not a large discrepancy in terms 
of magnitude. However, minor unknowns in terms of skin stresses 
for hot structures can have serious consequences. Skins are very 
susceptible to buckling because of prevailing compressive thermal 
stresses and stiffness degradation due to high temperature. The 
high temperature of the skins also creates the environment for 
creep problems if the stresses are not accurately known. The 
calculated thermal stresses in the skin areas exhibited distinct 
sensitivity to temperature smearing variations. It is felt that 
the problems with predicting skin thermal stresses with the 
NASTRAN models was more likely attributed to inadequate tempera- 
ture inputs rather than modeling insufficiencies. 

The ASTEEL model was developed with 97 bar and shear panel 
elements. The model contained approximately 285 elements per 
square meter (27 elements per square foot) of surface area. 
Approximately 147 temperature inputs per square meter (14 per 
square foor) of surface area were used to develop predicted 
thermal stresses. The CSTEEL model had 183 bar and shear panel 
elements. There were approximately 540 elements per square meter 
(50 elements per square foot) of surface area. The element 
density was considered adequate to predict the distribution of 
thermal stress and the time histories of thermal stress reliably. 
Approximately 224 temperature inputs per square meter (21 per 
square foot) of surface area were required to develop predicted 
thermal stresses. 


CONCLUDING REMARKS 


Several important problems were addressed in this paper. 

A structure fabricated of dissimilar materials was investigated 
in terms of correlating measured and calculated thermal stresses. 
Additional information was added to a growing data base from which 
estimates of finite element model complexities can be made with 
respect to thermal stress analysis. The manner in which tempera- 
tures were smeared to the finite element grid points was examined 
from the point of view of the impact on thermal stress calculations. 

The general comparison of calculated and measured thermal 
stresses is quite good and there is little doubt that the finite 
element approach provided by NASTRAN results in correct thermal 
stress calculations. Discrepancies did exist between measured 
and calculated values in the skin and the skin/frame junctures. 

The problems with predicting skin thermal stresses were attributed 


7 



to inadequate temperature inputs to the structural model rather 
than modeling insufficiencies. The discrepancies occurring at 
the skin/frame juncture were most likely due to insufficient 
modeling elements rather than temperature problems. 

In some areas, calculated thermal stresses varied consid- 
erably due to different temperature smearing approaches. It is 
apparent from the results of this paper, that thermal stresses 
calculated using the NASTRAN models are very sensitive to minor 
variations in the temperature input to the grid points. 

The more sophisticated NASTRAN model (CSTEEL) was developed 
with 183 bar and shear panel elements. The model contained 
approximately 540 elements per square meter (50 elements per 
square foot) of surface area. The element density was considered 
adequate to predict the distribution of thermal stress and the 
time histories of thermal stress reliably. Approximately 224 
temperature inputs per square meter (21 per square foot) of 
surface area were required to develop predicted thermal stresses. 


Dryden Flight Research Center 

National Aeronautics and Space Administration 
Edwards, Calif., June 6, 1979 
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Figure 6. Time history of altitude and Mach number for a hypersonic flight. 
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Figure 7. Test structure showing location of structural model. 
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Figure 8. Test structure showing location of structural model. 
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Figure 9. Grid point numbering system for NASTRAN models. 
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Figure 10. Measured temperatures and comparison of measured and 
NASTRAN -calculated thermal stresses for hypersonic heating simulation. 
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Figure 10. Continued. 
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Figure 10. Continued . 
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Figure 10. Concluded. 
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Figure 11. Comparison of measured and NASTRAN -calculated thermal stress distribution for the hyper- 
sonic profile. 
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Figure 12. Location of thermocouples with respect to the NASTRAN structural 
model . 
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Figure 13. Temperature at each of the skin heating control zones 
for time = 4 minutes. 
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Figure 14. Illustration of longitudinal temperature smearing. 
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Figure 14. Concluded . 
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Figure 15. Illustration of the lateral temperature smearing. 
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elements rather than temperature problems. In some areas, cal- 
culated thermal stresses varied considerably due to different 
temperature smearing approaches. It is apparent from the results 
of this paper, that thermal stresses calculated using the NASTRAN 
models are very sensitive to minor variations in the temperature 
input to the grid points. 
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